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VerticEil  shear  plus  horizontal  stretching  as  a route 
to  mixing 

Peter  H.  Haynes 

Department  of  Applied  Mathematics  eind  Theoretical  Physics,  University  of  Cambridge,  UK 


Abstract.  The  combined  effect  of  vertical  shear  and  horizontal  stretching 
leads  to  thin,  sloping  structures  in  tracer  fields,  whose  vertical  length  scale 
is  much  smaller  than  their  horizontal  length  scale.  These  structures  are 
then  vulnerable  to  vertical  mixing  processes.  This  effect  needs  to  be  taken 
into  accoimt  when  interpreting  the  horizontal  structure  of  oceanic  (and 
atmospheric)  tracers  and  may  explain  recent  oceanographic  observations. 


1.  Introduction 

In  both  atmosphere  and  ocean  it  is  useful  to  divide 
velocity  fields  into  isentropic  or  isopycnal  parts,  which 
move  fluid  parcels  along  potential  temperature  or  po- 
tential density  smfaces,  and  diabatic  or  diapycnal  parts, 
which  move  fluid  parcels  across  such  surfaces.  This  is 
because  diabatic  or  diapycnal  motion  can  be  accom- 
plished only  by  non-conservative  physical  processes  such 
as  radiative  heating  (in  the  atmosphere)  and  molecular 
diffusion  of  temperature  enhanced  by  three-dimensional 
turbulence  (in  atmosphere  and  ocean).  In  large  parts  of 
the  atmosphere  and  ocean,  such  processes  are  weak  in 
some  average  sense,  so  that  diabatic  or  diapycnal  veloc- 
ities must  be  much  smaller  than  isentropic  or  isopycnal 
velocities. 

The  process  of  stirring  may  therefore  be  regarded  as 
layerwise  two-dimensional,  in  the  sense  that  differential 
advection  by  the  quasi- horizontal  isentropic/isopycnal 
velocity  field  acts  independently  on  each  isentropic  or 
isopycnal  surface  to  distort  tracer  fields  into  complex 
geometrical  configurations.  If  atmospheric  and  oceanic 
flows  were  purely  two-dimensional,  similar  to  many  nu- 
merical simulations  or  laboratory  experiments,  then 
stirring  would  simply  strengthen  horizontal  spatial  gra- 
dients of  tracer  fields  imtil  horizontal  diffusion  became 
competitive  with  advection  and  mixing  occured.  But 
in  layerwise  two-dimensional  atmospheric  and  oceanic 
flows,  horizontal  advection  that  varies  in  the  vertical 
may  also  act  to  increase  vertical  gradients,  thereby  en- 
hancing the  effects  of  vertical  diffusion,  or  more  compli- 
cated vertical  mixing  processes.  The  fact  that  vertical 
scales  are  almost  always  observed  to  be  considerably 
smaller  than  horizontal  scales  suggests  that  mixing  dif- 
fusion rather  than  horizontal  mixing  is  the  dominant 
process. 

Most  previous  theoretical  investigations  of  stirring 


and  mixing  of  tracers  have  focussed  on  models  of  quasi- 
isotropic flows,  either  in  two  or  three  dimensions.  The 
considerations  above  suggest  that  the  most  relevant 
model  for  the  atmosphere  and  the  ocean  is  an  anisotropic 
model  in  which  the  vertical  velocity  is  zero,  but  the 
horizontal  velocities  depend  on  the  vertical  coordinate. 
Haynes  and  Anglade  (1997,  hereafter  HA97)  used  such 
an  anisotropic  model  to  study  the  enhancement  of  ver- 
tical gradients  by  horizontal  stirring  plus  vertical  shear. 
This  is  a relatively  simple  extension  to  previous  work, 
but  essential  if  theoretical  results  are  to  be  relevant  to 
real  flows.  The  HA97  work  was  motivated  primarily  by 
the  stratosphere.  This  article  reviews  and  expands  on 
some  of  the  eirguments  in  that  paper,  focussing  in  par- 
ticular on  the  implications  for  stirring  and  mixing  in 
the  ocean. 

The  model  to  be  presented  is  based  on  the  paradigm 
of  ‘chaotic  advection’,  that  flows  that  are  smoothly 
varying  in  space  and  time  may  lead  to  patterns  of  ad- 
vected  tracer  that  are  spatially  highly  complex.  One 
of  the  requirements  on  a chaotic  advection  flow  is  that 
there  is  a useful  division  of  scales  between  the  smallest 
active  scale  in  the  advecting  velocity  field  and  the  dif- 
fusive scede  (or  more  generally  scale  at  which  genuine 
mixing  processes  act).  The  case  of  three-dimensional 
turbulence,  for  example,  does  not  fall  into  this  cate- 
gory, since  the  velocity  field  has  complex  spatial  struc- 
ture at  scales  down  to  the  diffusive  scale  and  it  may 
be  argued  that  the  stirring  of  tracer  features  at  any 
scale  larger  than  the  diffusive  scale  is  dominated  by 
the  flow  at  the  same  scale  (rather  than  at  the  large 
scale).  In  fact  the  chaotic  advection  paradigm  is  rele- 
vEint  in  three-dimensionally  turbulent  flows  for  tracers 
with  diffusivity  much  less  than  the  momentum  diffusiv- 
ity  (large  Schmidt  number).  The  tracer  field  may  then 
have  non-trivial  structure  on  a spatial  scale  smaller  than 
the  Kolmogorov  scale,  with  the  flow  at  the  Kolmogorov 
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scale  providing  the  ‘lairge-scale’  stirring.  This  is  often 
called  the  Batchelor  regime  of  turbulence  and  has  been 
studied  as  an  important  special  case  in  recent  theoret- 
ical work  on  tracer  fields  in  turbulent  flows  (e.g.  see 
Shraiman  and  Siggia,  2000  and  references  therein). 

In  the  stratosphere  it  seems  to  be  useful  to  consider 
tracer  fields  to  be  determined  by  the  stirring  effect  of  the 
large-scale  flow  (varying  on  spatial  scales  of  hundreds 
to  thousands  of  kilometers)  and  by  the  mixing  effect  of 
localised  patches  of  three-dimensional  turbulence,  with 
vertical  scales  of  a hundred  meters  or  so  and  horizontal 
scales  of  a few  kilometers.  There  is  no  strong  evidence 
of  active  stirring  by  eddies  on  scales  from  a few  hundred 
down  to  tens  of  kilometers,  and  in  this  range  of  scales 
the  chaotic  advection  paradigm  is  therefore  relevant. 

Similarly  in  the  ocean  one  might  argue  that  mesoscale 
eddies  with  scales  of  tens  of  kilometers  are  the  dominant 
part  of  the  flow  in  stirring  the  tracer  field  at  scales  from 
a few  tens  of  kilometers  down  to  the  scale  of  small-scale 
mixing  events. 

The  structure  of  the  paper  is  as  follows.  Section  2 de- 
scribes a simple  model  problem  of  steady  flow  including 
both  horizontal  stretching  and  vertical  shear,  focussing 
on  the  spreading  of  a tracer  from  a point  release  in  such 
a flow.  Section  3 discusses  some  aspects  of  tracer  evolu- 
tion in  more  general  random  flows.  Section  4 considers 
application  of  these  ideas  to  the  observations  of  Ledwell 
et  al.  (1993)  (hereafter  L93).  Section  5 gives  a brief 
discussion  of  the  possible  competing  role  of  mixing  by 
double-diffusive  intrusions.  Section  6 describes  some 
possible  future  lines  of  research. 

2.  Steady  flow  model 

The  simplification  allowed  by  the  chaotic  advection 
paradigm  is  that  there  is  a finite  scale  below  which  the 
flow  may  be  considered  as  a linear  function  of  the  space 
coordinates,  i.e.  the  velocity  field  may  be  expanded 
as  a (multidimensioneil)  Taylor  series  about  some  ref- 
erence location  and  only  the  linear  term  in  the  series 
retained.  The  effects  of  small-scale  mixing  events  may 
represented  by  a diffusion  term  (but  more  sophisticated 
representations  are  also  possible). 

We  follow  HA97  in  considering  first  a steady  linear 
flow  that  has  the  two  ingredients  of  horizontal  strain,  F, 
plus  vertical  shear,  A.  In  much  of  the  atmosphere  and 
ocean  it  is  relevant  to  consider  the  parameter  regime 
A » r,  i.e.  vertical  shear  much  larger  than  horizontal 
strain.  This  may  be  argued  from  the  observed  fact  that 
horizontal  length  scales  tend  to  be  larger  them  verti- 
cal scales  or  by  appealing  to  theoretical  arguments  that 
the  ratio  of  horizontal  to  vertical  scales  is  of  order  f/N 
(Prandtl’s  ratio),  where  / is  the  Coriolis  parameter  and 


N the  buoyancy  frequency.  f/N  is  0(100)  in  the  at- 
mosphere and  in  the  oceanic  thermocline  and  perhaps 
reduces  to  0(10)  in  the  deep  ocean. 

The  flow  is  taken  to  have  components  (Fx,  —Ty  -|- 
Az,  0)  in  the  x,  y and  z directions  respectively,  where 
(x,  y,  z)  are  Cartesian  coordinates,  with  z vertical.  The 
steadiness  assumption  implies  that  the  horizontal  strain 
F and  the  vertical  shear  A are  constants.  The  equation 
for  the  evolution  in  this  flow  of  a tracer  with  concentra- 
tion x(a:,  y,  2,  t)  is  therefore 


^+rx^  + (Az-Fu)^  = K(^-b^-f^)  (1) 
dt^  dx^^  ^’dy  dy^^  dz^’^ 


where  k is  the  diffusivity,  assiuned  constant. 

HA97  considered  sinusoidal  solutions  of  this  equa- 
tion. Another  useful  solution  of  this  equation,  and  in- 
deed of  the  equation  for  tracer  concentration  in  any  flow 
where  velocity  components  are  linear  functions  of  space 
coordinates  (including  time-dependent  flows),  is  an  el- 
lipsoidal Gaussian  solution,  i.e.  an  exponential  func- 
tion of  a negative  definite  quadratic  function  oi  x,  y 
and  2.  This  is  conveniently  obtained  by  considering  the 
second-order  spatial  moments  of  the  solution,  such  as 
mil  = Jx^x^/  / The  six  moments 

are  sufficient  to  determine  the  quadratic  function  and 
the  constraint  that  the  total  amount  of  tracer  remains 
constant  then  determines  the  coefficient  of  the  exponen- 
tial function.  In  the  above  flow,  the  equations  for  the 
moments  may  be  shown  to  be 


1 

dt 


d 

at 

(2) 

— A,771xz 

dt  ^ 

(3) 

d 

~77lxz  — I^TT^xz 

(4) 

rriyy  = ~2Tmyy  -f  2Amyx  + 2/C 

(5) 

(6) 

d 

-m„  = K 

(7) 

These  equations  may  be  solved  straightforwardly. 
For  brevity  we  consider  the  case  where  all  moments 
tend  to  zero  as  t -)•  O'*',  corresponding  to  a point  re- 
lease. Then  the  solutions  are 


mil  = =e 


2rt 


(8) 


"^xy  — ^ 


(9) 
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rr^yy  = ^ {2r*  + 4e-"‘  - - 3}+^  {l  - 

(10) 

icA 

^3/*  = {ri- l + e“^‘}  (11) 

rrizz  = 2nt  (12) 

First  consider  the  case  where  A = 0,  i.e.  the  case  of 
purely  two-dimensional  flow.  Thenm^a;  = K(e^^*  — 1)/F 
and  niyy  = k(1  — e~^'^*)/r.  For  small  times,  t F~^, 
diffusion  dominates  and  both  rrixx  and  niyy  increase  lin- 
early with  time.  When  < ~ F~^  advection  by  the  hor- 
izontal strain  flow  becomes  competitive  with  diffusion. 
The  tracer  elongates  in  the  direction  of  the  stretch- 
ing axis,  so  that  for  t » F"^,  rrixx  ~ whilst 

in  the  direction  of  the  compression  axis  a steady-state 
balance  between  strain  and  diffusion  is  achieved,  with 
rriyy  ~ k/F.  At  all  times  the  only  vertical  transport  is 
through  diffusion,  so  rrixx  increases  linearly  with  t at  all 
times. 

What  is  the  effect  of  adding  vertical  shear,  i.e.  tak- 
ing 0 (with  A ^ F)?  The  elongation  of  the  tracer 
along  the  stretching  axis  of  the  horizontal  strain  axis  is 
imchanged.  However,  the  spreading  of  the  tracer  in  the 
direction  of  the  compression  axis  is  rather  different.  Ex- 
amination of  the  expression  for  niyy  above  shows  that 
there  are  three  regimes.  For  t < A~^  (regime  I)  there  is 
purely  diffusive  spreading  as  above.  For  A“^  < F~^ 

(regime  II)  the  diffusive  spreading  of  the  tracer  in  the 
vertical  allows  the  vertical  shear  to  augment  the  spread- 
ing in  the  horizontal,  so  that  rriyy  ~ kA’^I^.  Finally 
when  t ~ F“^  horizontal  advection  begins  to  inhibit 
the  horizontal  spreading,  so  that  for  f>F~^  (regime 
III),  rriyy  f\j  KtA^/F^.  These  three  regimes  are  depicted 
graphically  in  Figure  1.  Note  that  rriyy  continues  to  in- 
crease with  time,  in  contrast  to  the  case  with  A = 0. 
In  regime  III  rriyy  ~ Arriyx/T  ~ A'^rrixz/T^  suggest- 
ing that  the  y and  z scales  of  the  tracer  patch  are  in 
the  ratio  A/F.  Indeed  what  happens  is  that  the  tracer 
patch  forms  a sheet  that  slopes  at  angle  F/A  (from 
the  horizontal)  in  the  (y,  z)  plane.  The  extent  of  the 
sheet  in  the  vertical  direction  is  («;f)^/^,  therefore  the 
extent  of  the  sheet  in  the  horizontal  (y)  direction  is 
A(/ct)^/^/F.  This  explains  the  dominant  behaviour  of 
rriyy,  rriyx  and  rrixz-  In  order  to  deduce  the  spread  of 
the  tracer  on  a given  horizontal  surface  it  is  useful  to 
consider  rriyy  — 2Arriyx/T  + rrixzA!^ , i-e.  the  moment 
of  (y  — Az/F)^.  It  is  straightforward  to  show  that 

rriyy  - 2Amyz/T  -t-  rrizzA^ = 

K A^ . . ^ ft  —9vt\ 

(13) 


(Ko-a  (I) 

^ At  (Kt)'^^ 


(A/rr-(Kt)>« > 


Figure  1.  The  spread  of  tracer  in  the  (y,z)  in  each  of 
regimes  I {t<  A"^),  regime  II  (A"^<  t<  F~^)  and  regime 

Iii(r-i<t). 

Thus  the  spread  of  the  tracer  on  a given  horizon- 
tal surface  is  (k/F)^/^A/F.  The  effect  of  the  vertical 
shear  is  therefore  to  increase  the  equilibrium  width  of 
the  tracer  patch  by  a factor  A/F. 

In  summary,  the  effect  of  the  vertical  shear  is  to  lead 
to  sloping  structures  in  the  tracer  field  with  aspect  ratio 
(ratio  of  horizontal  length  scale  to  vertical  length  scale) 
a = A/F.  The  equlibrium  width  of  these  structures 
is  larger  by  a factor  a than  would  be  the  case  without 
vertical  shear.  One  might  say  that  it  is  as  if  a horizontal 
diffusivity  na^  were  acting. 

3.  Time-dependent  flow  models 

In  a flow  that  varies  in  space  and  time,  the  velocity 
gradient  encountered  by  a fluid  parcel  changes  with  time 
as  the  particle  moves  through  the  flow.  The  evolution  of 
small-scale  tracer  features  is  therefore  governed  by  the 
equation  for  tracer  evolution  in  a time-dependent  linear 
flow.  In  the  study  of  turbulence  there  is  a long  tradition 
of  considering  the  effect  of  a linear  flow  that  varies  ran- 
domly in  time  (using  so-called  ‘random-straining’  mod- 
els). This  approach  has  been  used  most  recently  by 
various  authors  to  give  a rather  complete  description  of 
the  statistics  of  tracer  fields  in  Batchelor-regime  turbu- 
lence. [See  Balkovsky  and  Fovxon,  1999  and  Falkovich 
et  al.,  2001  for  further  details  of  this  work.] 

The  random-straining  model  most  relevant  to  the  at- 
mosphere and  ocean  is  non-standard  in  that  the  statis- 
tics of  the  velocity  gradient  tensor  is  not  isotropic  (since 
the  vertical  velocity  is  zero,  but  the  vertical  gradient  of 
the  horizontal  velocity  is  not).  HA97  considered  such 
random-straining  models  and  showed  that  the  impor- 
tant predictions  of  the  steady  flow  model  presented  in 
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Section  2 carried  over,  in  the  sense  that  the  tracer  field 
was  predicted  to  form  sloping  sheets,  with  an  aspect  ra- 
tio a that  depended  on  the  statistics  of  the  horizontal 
strain  and  vertical  shear  fields. 

We  review  some  of  the  HA97  results  below,  empha- 
sising the  dependence  of  a on  the  different  flow  parame- 
ters. One  particular  goal  is  to  identify  conditions  under 
which  a is  not  simply  equal  to  the  aspect  ratio  A/T  of 
the  flow  itself. 

It  is  convenient  to  consider  the  normalised  gradient  of 
the  tracer  field,  defined  by  k(t)  = evaluated  at 

the  position  x = X(t)  following  a fluid  parcel,  k might 
be  considered  as  a local  wavenumber  of  the  tracer  field 
in  the  case  where  the  tracer  field  varies  rapidly  in  space 
compared  to  the  velocity  field.  It  may  be  shown  from 
the  gradient  of  the  tracer  advection  equation  that 

^ = -(Vu).k  (14) 

where  the  scalar  product  on  the  right-hand  side  applies 
to  the  second  index  of  the  tensor  Vu.  The  evolution  of  k 
therefore  depends  on  the  time  series  of  Vu  (encountered 
following  the  point  X(t)).  In  a random-straining  model 
the  effect  of  turbulent  flow  on  gradients  of  tracer  (or 
line  elements)  is  considered  as  a kinematical  problem, 
in  which  the  tensor  Vu  is  taken  to  be  given  by  a random 
time  series,  without  addressing  questions  of  the  flow 
dynamics. 

For  the  layerwise  two-dimensional  flows  of  interest 
here,  the  matrix  W of  components  of  the  velocity  gra- 
dient tensor,  with  Wij  = dui/dxj,  has  W31  = W32  = 
W33  = 0.  (The  index  3 is  taken  to  correspond  to 
the  z coordinate.)  It  is  useful  to  write  k = (k^^^m), 
with  k^^^  the  horizontal  wavenumber  and  m the  verti- 
cal wavenumber,  and  define  to  be  the  2nd  rank 
tensor  with  components  Wn,  W21,  W12  and  W22,  and  A 
to  be  the  vector  with  components  (W13,  W23).  It  then 
follows  that 

^ (15) 

and 

5 (16) 

In  what  follows  the  summation  convention  is  used  for 
repeated  suffices  and  the  suffices  run  through  the  val- 
ues {1,2}.  Equation  (15)  is  just  that  for  the  tracer 
wavenumber  in  twodimensional  (horizontal)  flow.  Note 
that  it  may  be  solved  independently  of  (16).  Thus  vari- 
ation in  the  z-direction  makes  no  difference  to  the  evolu- 
tion of  the  horizontal  wavenumber  vector.  However  (16) 
shows  that  the  vertical  wavenumber  evolves  through  the 
vertical  shear  (in  the  horizontal  flow)  acting  on  the  hor- 
izontal wavenumber  vector. 


To  formulate  a suitable  random-straining  model  we 
assume  that  each  of  the  six  non-zero  components  of  the 
matrix  W is  a realisation  of  a random  function  of  time, 
with  average  (over  all  realisations)  zero.  We  also  assume 
that  the  wavenumber  at  t = 0 is  randomly  chosen  and 
is  statistically  independent  of  the  Wij.  Equation  (16) 
may  be  integrated  to  give 

m(f)  = m(0)-h  f Aj{t')kf\t')dt'  (17) 

Jo 

and  then  squaring,  and  taking  the  ensemble  average,  it 
follows  that 

(m(t)2)  = (m(0)^)+ 

/‘df'  rdt"(A,(t')A:f(t')Afe(f")^r^(<")>,  (18) 

Jo  Jo 

where  (.)  denotes  the  ensemble  average,  over  all  reali- 
sations of  the  flow  field  and  all  realisations  of  the  initial 
wavenumber  vector. 

Useful  quantitative  estimates  are  possible  if  a num- 
ber of  further  assumptions  are  made  concerning  the 
statistics  of  the  velocity  gradient  tensor.  The  first  is 
that  the  components  of  A are  independent  of  all  com- 
ponents of  i.e.  that  the  vertical  shear  is  statis- 

tically independent  of  the  horizontal  deformation.  It 
follows  that  A is  also  independent  of  k^^^  so  that  av- 
erages of  terms  involving  A and  terms  involving  k^^^ 
may  be  taken  separately  in  (18).  The  second  impor- 
tant assumption  is  stationarity,  from  which  it  follows 
||(Aj(f')Afc(t"))||  ~ A^g{aA\t'  — t"\),  for  some  function 
g,  where  a a is  an  inverse  correlation  time  for  the  vertical 
shear  and  A^  = (A^). 

In  order  to  estimate  the  integral  appearing  in  (18) 
it  is  also  necessary  to  have  information  on  the  evolu- 
tion of  the  horizontal  wavenumber.  This  follows  ex- 
actly as  in  the  random  straining  models  of  isotropic 
twodimensional  flow  considered  by  Kraichnan  (1974) 
and  others.  The  basic  predictions  of  these  theories  are 
that  the  horizontal  wavenumber  increases  exponentially 
in  time,  at  a rate,  S say,  governed  by  the  statistics 
of  the  horizontal  strain  field  and  depending  in  partic- 
ular on  the  root  mean  square  rate  of  strain,  T say, 
and  on  the  inverse  correlation  time  ah  for  the  hori- 
zontal str2iin  field.  Detailed  analysis  of  explicit  models 
[Kraichnan,  1974;  Chertkov  et  al.  , 1995;  HA97)  sug- 
gests that  5 ~ rmin(r/cT/i,  1). 

Again,  making  suitable  stationarity  assumptions,  this 
suggests  that  (<”))! I ~ kQe^^^''^*"^h[ij.\t'  - 

t"|),  where  kg  = (k^^)(0)^).  The  inverse  time  scale  /j, 
and  the  precise  form  of  the  function  h(.),  depend  on 
the  statistics  of  the  horizontal  strain  field,  but  explicit 
calculation  in  various  models  suggests  that  /it  ~ 5. 
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Substituting  these  estimates  into  (18)  it  follows  that, 
for  large  times,  when  the  second  term  on  the  right-hand 
side  dominates  the  first. 


imn{l,  — }. 


(19) 


Note  that  the  dominant  contribution  to  the  first  inte- 
gral, over  t",  performed  in  (18)  comes  from  a neighbour- 
hood of  t"  = t of  size  min{5“^,a)(^^}  and  that  to  the 
second  integral,  over  t',  comes  from  a neighbourhood  of 
t’  = t ol  size  S~^. 


The  ratio  a of  horizontal  length  scale  to  vertical 
length  scale,  or,  equivalently,  the  ratio  of  vertical  wavenum- 
ber to  horizontal  wavenumber  may  therefore  be  esti- 
mated as 


--  ^max(l,  )H20) 

i 1 CT\  Oh 


with  the  last  estimate  following  from  the  estimate  for 
S. 

This  semi-quantitative  analysis  suggests  that  the  ver- 
tical wavenumber  increases  exponentially  at  the  same 
rate  as  the  horizontal  wavenumber,  and  that  the  ratio  a 
is  equal  to  the  ratio  of  vertical  shear  to  horizontal  strain, 
A/r,  multiplied  by  a number  that  depends  on  Oh/T  and 
cta/T.  These  predictions  agree  with  explicit  calculations 
based  on  suitably  formulated  random-straining  models. 
(See  HA97  for  more  details.) 

Figure  2 summarises  the  variation  of  the  aspect  ra- 
tio a with  the  parameters  Oh  and  a a..  Note  that  a is 
anomalously  large  (i.e.  greater  than  A/F)  when  the  hor- 
izontal strain  varies  on  a timescale  that  is  shorter  than 
both  the  inverse  horizontal  strain  rate  and  the  timescale 
of  variation  of  the  vertical  shear.  On  the  other  hand, 
a is  anomalously  small  (i.e.  less  than  A/F),  when  the 
vertical  shear  varies  on  a timescale  that  is  shorter  than 
both  the  inverse  horizontal  strain  rate  and  the  timescale 
of  variation  of  the  horizontal  strain.  As  might  be  ex- 
pected, in  the  limit  where  cr/i  ■C  F and  cta  ^ F the 
estimate  for  a agrees  with  that  obtained  for  the  steady 
flow  considered  in  the  previous  section. 

Note  that  in  this  section  we  have  not  considered  the 
effects  of  diffusion  exphcitly.  Nonetheless,  it  is  plausible 
that  diffusion  acts  in  a similar  way  to  that  deduced 
from  the  steady  flow  model  in  Section  2,  i.e.  diffusion 
K acts  on  structures  in  the  tracer  field  that  slope  with 
eispect  ratio  a and  the  result,  e.g.  in  achieving  a balance 
between  diffusion  and  horizontal  straining,  is  that  it  is 
as  if  there  is  a horizontal  diffusivity  with  magnitude 
Ka^.  This  is  supported  (with  certain  hmitations)  by 
explicit  calculations  in  Vanneste  and  Haynes  (2001). 


Figure  2.  Scaling  of  (where  a.  is  aspect  ratio  of  hori- 
zontal to  verticeJ  scale)  with  parameters  F (horizontal  strain 
rate),  A (vertical  shear),  Oh  (inverse  correlation  time  for  hor- 
izontal strain)  2ind  o\  (inverse  correlation  time  for  vertical 
shear).  The  thin  dashed  lines  delimit  the  regions  where  the 
different  estimates  hold.  The  thick  dashed  lines  delimit  the 
regions  where  a scales  as  the  2ispect  ratio  of  the  velocity  field 
(A/r)  (lower  left-hand  region),  is  larger  than  that  aspect  ra- 
tio. (right-hand  region)  and  is  smaller  them  that  aspect  ratio 
(upper  region). 


4.  Application  to  the  Ledwell  et  al. 
(1993)  observations 

In  the  tracer  release  experiment  reported  by  L93  the 
dispersion  of  a tracer  in  the  ocean  thermocline  was 
followed  over  a period  of  several  months.  The  verti- 
cal (cross-isopycnal)  diffusion  could  be  estimated  di- 
rectly and  a value  of  vertical  diffusivity  of  about 
lO^^m^s'^was  inferred.  In  the  horizontal  the  tracer 
was,  in  the  later  stages  of  the  experiment,  observed 
to  be  confined  to  thin  streaks,  whose  width  apparently 
reached  an  equilibrium  value  of  about  3 km.  This  was 
interpreted  as  the  stretching  out  of  the  tracer  patch  by 
the  mesoscale  eddy  field,  resulting  in  filaments  of  tracer 
whose  width  was  the  equilibrium  value  determined  by  a 
balance  between  horizontal  stretching  (which  tends  to 
reduce  the  width  of  the  filament)  and  horizontal  disper- 
sive effects,  perhaps  associated  with  small-scale  mixing 
processes  (which  tend  to  increase  the  width  of  the  fila- 
ment). Prom  the  observed  length  of  the  tracer  filament 
it  was  estimated  that  the  average  stretching  rate  expe- 
rienced was  about  3 x 10“’^s"^.  It  followed  that  if  the 
horizonted  dispersive  effects  could  be  represented  by  an 
effective  horizontal  diffusivity  ka,  it  must  have  a value 
of  about  3 m^s~^(ci  (3  km)^  x 3 x 10“^s“^).  Thus 
Kh/Ky  ~ 3 X 10®.  Previous  work  by  Young  et  al.  (1982), 
taking  account  of  the  combined  effect  of  vertical  mixing 
and  horizontal  advection  by  inertio-gravity  waves  had 
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suggested  kh/kv  ^ (N/f)^  ~ 1500,  i.e.  200  times  too 
small  to  account  for  the  L93  observations. 

HA97  noted  that  the  combined  effect  of  horizontal 
stretching  and  vertical  shear  could  account  for  the  value 
of  Kh/Ky  if  the  aspect  ratio  a was  about  500.  According 
to  the  analysis  presented  in  Section  3 this  requires  A ~ 
1.5  X 10~^s~^,  based  on  the  estimate  a ~ A/S  (which 
holds  if  ca  ^ S,  i.e.  the  correlation  time  for  the  vertical 
shear  must  be  larger  than  S~^).  This  value  of  vertical 
shear  (equivalent  to  15  cm  s"^  per  km)  does  not  seem 
indefensible  as  reahstic. 

If  the  models  described  in  this  article  are  to  be  rel- 
evant to  the  L93  observations,  one  might  ask  further 
whether  the  observed  morphology  of  the  tracer  is  consis- 
tent with  the  model  predictions.  First  it  is  clear  that  in 
a flow  with  flnite  length  scales  the  exponential  stretch- 
ing in  a single  direction  predicted  by  the  model  in  Sec- 
tion 2 must  break  down  when  the  largest  length  scale  of 
the  tracer  patch  (or,  at  this  stage,  tracer  filaments)  be- 
comes comparable  to  the  length  scale  of  the  flow.  What 
happens  is  simply  that  the  filaments  meander  on  the 
length  scale  of  the  flow  (as  observed  in  countless  labora- 
tory and  numerical  experiments).  Perhaps  more  crucial 
is  the  shape  of  tracer  filaments  in  cross-section.  The 
model  of  Section  2 predicts  that  this  will  ultimately  be 
highly  elongated,  so  that  the  tracer  is  actually  concen- 
trated within  sloping  sheets.  However,  this  applies  on 
time  scales  much  greater  than  T“^,  where  T is  the  strain 
rate  and  therefore,  in  the  steady  model,  the  stretching 
rate.  For  times  comparable  to  F“^  (see  Figure  1),  the 
tracer  distribution  is  not  sheet-like,  but  it  is  simply  the 
case  that  the  filament  cross  section  has  much  larger  hor- 
izontal extent  than  vertical  extent  (by  a factor  q).  The 
observing  period  reported  by  L93  is  such  that  Ft  ;$  6, 
so  that  sheet-like  features  are  not  necessarily  expected. 
[Note  that  the  extent  of  the  sheets  is  proportional  to 
(rt)i/2.] 

5.  Possible  effects  of  double  diffusion 

Garrett  (1982)  (hereafter  G82)  noted  that  the  stir- 
ring by  mesoscale  eddies  of  temperature  and  salinity 
along  isopycnal  surfaces  would  lead  to  large  local  gra- 
dients in  these  quantities  manifested  as  thermohaline 
fronts.  These  in  turn  might  lead  to  double-diffusive  in- 
trusions. G82  argued  that  the  sharpening  of  gradients 
of  temperature  and  salinity  by  stirring  would  be  halted 
by  the  mixing  effects  of  the  intrusions  when  the  growth 
rate  of  the  intrusions  was  equal  to  the  convergence  act- 
ing on  the  fronts,  i.e.  to  the  stretching  rate  associated 
with  the  stirring  process.  He  gave  the  formula  (based 
on  previous  theoretical  and  experimental  work) 

= 0.075^  (21) 


for  the  maximum  growth  rate  Xmax  of  the  intru- 
sions, where  g is  the  gravitational  acceleration,  p = 
p~^{dp/dS)r  is  the  rate  of  change,  at  constant  tem- 
perature, of  density  with  saHnity,  and  Sx  is  the  salinity 
gradient  in  the  neighbourhood  of  the  front.  He  argued 
that  if  the  front  had  width  Lf  then  the  approximate 
value  of  Sx  would  be  Sx  — SxLeddy/Lf  where  Sx  was 
the  large-scale  salinity  gradient  and  Leddy  was  the  typ- 
ical size  of  the  mesoscale  eddies.  If  s is  the  stretching 
rate  then  it  follows  that  Lf  is  given  by  the  formula 

Lf  = 0.075^^^^.  (22) 

Ns 

For  the  values  of  s and  Lf  observed  by  L93  (taking 
Lf  to  be  the  width  of  the  filaments)  and  assuming 
Leddy  ~ 100  km,  this  would  require  gPSx/N  2;  10“^s“^. 
This  is,  in  fact,  considerably  weaker  than  the  values  sug- 
gested by  G82  as  examples.  This  therefore  appears  to 
leave  open  the  possibihty  that  it  is  the  mixing  effects, 
along  isopycnal  surfaces,  of  double  diffusive  intrusions, 
rather  than  any  process  directly  related  to  ‘background’ 
vertical  mixing,  that  determines  the  width  of  the  fila- 
ments observed  by  L93.  On  the  other  hand  it  is  not 
at  all  clear  that  the  filaments  of  a tracer  injected  at  an 
arbitrary  location  will  match  the  regions  of  enhanced 
gradients  in  temperature  and  salinity  whose  distribu- 
tion is  set  on  the  large  scale  and  therefore  feel  the  same 
mixing  effects. 

6.  Discussion 

The  work  reported  in  HA97  may  be  extended  in  var- 
ious ways,  some  of  which  are  relevant  to  oceanographic 
considerations.  For  example,  Vanneste  and  Haynes 
(2001)  have  considered  the  effect  of  vertical  shear  on 
the  horizontal  wavenumber  spectrum  of  passive  tracers, 
in  particular  on  the  range  of  scales  where  the  effects  of 
diffusive  mixing  become  important.  In  particular  this 
work  highlights  the  hmitations  of  estimating  the  effec- 
tive horizontal  diffusivity  as  kq^.  The  recent  theoreti- 
cad  work  on  Batchelor-regime  turbulence  has  made  ex- 
plicit predictions  about  probability  density  functions  for 
tracer  concentrations,  tracer  concentration  differences 
(over  a finite  distance)  and  tracer  concentration  gradi- 
ents and  this  work  needs  to  be  extended  to  the  layerwise 
two-dimensional  case  if  it  is  to  be  applicable  to  real  at- 
mospheric and  oceanic  flows. 

Much  insight  into  atmospheric  tracer  distributions 
has  been  obtained  by  using  models  driven  by  observed 
velocity  fields,  or  velocity  fields  extracted  from  quasi- 
realistic  global  models.  Approaches  have  included  full 
numerical  integrations  of  the  tracer  evolution  equation 
(e.g.  on  an  isentropic  surface),  non-diffusive  reconstruc- 
tions of  tracer  fields  by  following  back  trajectories  to  an 
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initial  condition  and  calciilations  of  the  statistical  prop- 
erties of  large  numbers  of  trajectories,  e.g.  to  give  the 
distribution  of  finite-time  stretching  rates.  See,  for  ex- 
ample, the  papers  by  Schoeberl  and  Newman  (1995), 
Ngan  and  Shepherd  (1999)  and  Hu  and  Pierrehumbert 
(2001).  All  of  these  methods  have  substantial  savings 
over  integrations  of  a full  model  including  tracer  equa- 
tions and  d3mamical  equations  and  some  may  be  worth 
applying  in  the  oceanic  context.  Velocity  fields  might 
be  extracted  from  eddy-resolving  models  for  this  pur- 
pose, but  as  a first  step  it  might  be  worth  considering 
‘synthetic’  velocity  fields  that  are  generated  quite  artifi- 
cially to  have  a plausibly  reahstic  spatial  and  temporal 
structure  in  space  and  time. 
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